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ABSTRACT 

We present broadband observations of the afterglow and environment of the short GRB 1 1 1020A. An extensive 
X-ray light curve from Swift /XRT, XMM-Newton, and Chandra, spanning ~ 1 00 s to 10 days after the burst, reveals 
a significant break at St rs 2 days with pre- and post-break decline rates of a x . i — —0.78 and otx , 2 < —1-7, 
respectively. Interpreted as a jet break, we infer a collimated outflow with an opening angle of Oj ~ 3°-8°. The 
resulting beaming-corrected y-ray (10-1000 keV band) and blast-wave kinetic energies are (2-3) x 10 48 erg and 
(0.3-2) x 10 49 erg, respectively, with the range depending on the unknown redshift of the burst. We report a radio 
afterglow limit of <39 //Jy (3a) from Expanded Very Large Array observations that, along with our finding that 
fc < Vx, constrains the circumburst density to no ~ 0.01-0. 1 cm -3 . Optical observations provide an afterglow limit 
of i > 24.4 mag at 18 hr after the burst and reveal a potential host galaxy with i ~ 24.3 mag. The subarcsecond 
localization from Chandra provides a precise offset of O'.' 80 ± O'.' 1 1 (la ) from this galaxy corresponding to an offset 
of 5-7 kpc for z — 0.5-1. 5. We find a high excess neutral hydrogen column density of (7.5 ± 2.0) x 10 21 cm -2 
(z = 0). Our observations demonstrate that a growing fraction of short gamma-ray bursts (GRBs) are collimated, 
which may lead to a true event rate of >100-1000 Gpc” 3 yr 1 , in good agreement with the NS-NS merger rate 
of ~200-3000 Gpc 3 yr” 1 . This consistency is promising for coincident short GRB-gravitational wave searches in 
the forthcoming era of Advanced LIGO/ VIRGO. 
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1. INTRODUCTION 

Observations of the temporal and spectral evolution of short- 

duration gamma-ray burst (GRB, Tg 0 < 2 s; Kouveliotou et al. 

1993) afterglows are crucial to our understanding of the basic 
properties of these events: their energetics, parsec-scale envi- 
ronments, and geometries. From observations over the past 7 
years, we now know that short GRBs have isotropic-equivalent 
energies of ~10 5(, -10 52 erg (Berger 2007) and circumburst den- 
sities of — 10 fi to 1 cm” 3 (Soderberg et al. 2006; Panaitescu 
2006; Stratta et al. 2007; Perley et al. 2009b; Berger 2010; 
Fong et al. 2011); however, these ranges are based on only 
a handful of events. The geometry, or degree of collimation, 
is the least constrained property but is of particular interest 
because it directly affects the true energy scale and event 
rates. These parameters aid our understanding of the explo- 
sion physics, the nature of the progenitors, and the potential 
detectability of short GRBs as gravitational wave sources. In 
particular, knowledge of the true energy scale may constrain 
the mechanism of energy extraction from the central engine and 
the ejecta composition: vv annihilation powering a baryonic jet 
(Jaroszynski 1993; Mochkovitch et al. 1993) or magnetohydro- 
dynamic (MHD) processes in a magnetically dominated outflow 
(Blandford & Znajek 1977; Rosswog et al. 2003). Significant 
improvement on the short GRB observed rate of > 1 0 Gpc” 3 yr” 1 
(Nakar & Granot 2007) will have a critical impact on estimates 
for coincident short GRB-gravitational wave detections in the 
era of Advanced LIGO/ VIRGO (Abadie et al. 2010). 

The opening angles ( 0j ) of GRBs can be inferred from 
temporal breaks in the afterglow light curves (“jet breaks”), 


which occur at the time, tj, when the Lorentz factor of the 
outflow is r(tj') ~ I /Oj ; a later break corresponds to a wider 
opening angle (Sari et al. 1999; Rhoads 1999). Jet breaks in the 
light curves of long-duration GRBs have led to an opening angle 
distribution with a range of — 2°-20° and a median of 7° , leading 
to beaming-corrected energies of E y = [1 — cos(0 ; j | /L / j s0 ~ 
10 5 °-10 51 erg (Bloom et al. 2003; Frail et al. 2001; Friedman & 
Bloom 2005; Kocevski & Butler 2008; Racusin et al. 2009). For 
short GRBs, there is mounting theoretical (Eichler et al. 1989; 
Narayan et al. 1992) and observational (Fong et al. 2010; Berger 
2010; Church et al. 2011) evidence that the progenitors are 
NS-NS/NS-BH mergers, and numerous simulations of post- 
merger black hole (BH) accretion have predicted collimated 
outflows with 0j ~ 5°-20° (Popham et al. 1999; Aloy et al. 
2005; Rosswog 2005; Rezzolla et al. 201 1) up to several tens of 
degrees (Ruffert & Janka 1999b; Popham et al. 1999; Rezzolla 
et al. 2011). 

However, the detection of jet breaks in the afterglow light 
curves of short GRBs has proved to be challenging. They can 
in principle be measured from optical or radio observations, but 
there are several caveats that make this particularly difficult for 
short GRBs. First, the brightness of optical and radio afterglows 
is sensitive to the circumburst densities, which are typically 
low, ~10 2 cm” 3 (Soderberg et al. 2006). Indeed, of nearly 70 
short bursts detected by Swift, only 2 radio afterglows have been 
detected over the past 7 years (Berger et al. 2005; Soderberg et al. 
2006; Chandra & Frail 2011). Similarly, only -30% of Swift 
bursts have detected optical afterglows, with a typical brightness 
at <1 day of »23 mag (Berger 2010; Fong et al. 2011), 
making long-term temporal monitoring nearly impossible with 
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ground-based facilities. Second, in the optical band there can 
be significant contamination from the host galaxies, which are 
generally brighter than the afterglows at >1 day (Berger 2010). 

On the other hand, the X-ray afterglow brightness is inde- 
pendent of the circumburst density (as long as the density is 
>10“ 5 cm~ 3 and hence v c > v x ; Granot & Sari 2002), and host 
contamination is not an issue. In addition, the well-sampled 
Swift [X RT light curves from ~1 minute to ~1 day provide an 
unambiguous baseline against which we can measure a subse- 
quent break. Therefore, it is no surprise that the X-rays enabled 
the discovery of the first jet break in a short GRB. The X-ray 
afterglow light curve of GRB 05 1221 A exhibited a break at 
~5 days, leading to Qj ~ 7° (Soderberg et al. 2006; Burrows 
et al. 2006). Similarly, Chandra observations of GRB 050724A 
out to 22 days placed a meaningful lower limit of 0j > 25° 
(Grape et al. 2006), consistent with a spherical explosion. Tem- 
poral breaks on timescales of < few hours were observed in 
the XRT light curves of GRBs 061201 (Stratta et al. 2007) and 
0905 10 5 (De Pasquale et al. 2010). If these are interpreted as jet 
breaks, they would lead to Qj ~ 1°; however, they also match 
the timescale and behavior of early breaks in long GRBs, which 
are not due to collimation (Nousek et al. 2006; Zhang et al. 
2006; Liang et al. 2007). Finally, there is tentative evidence 
for beaming in the light curves of GRBs 050709 (Fox et al. 
2005) and 061210 6 (Berger 2007); however, these are based on 
sparsely sampled light curves without a definitive break (e.g., 
Watson et al. 2006). No other jet breaks in the light curves of 
unambiguous short GRBs have been reported to date, 7 and the 
lack of jet breaks in Swift /XRT observations out to ~l-2 days 
can provide only weak lower bounds of 0j > 2°-6° (revised 
from Coward et al. 2012 with more realistic density values; see 
Section 4.2). 

Against this backdrop, we present here the discovery of a 
break in the X-ray light curve of GRB 1 1 1020A at St 2 days, 
best explained as a jet break. We also present contemporaneous 
radio and optical limits on the afterglow, allowing a character- 
ization of the broadband synchrotron spectrum and constraints 
on the energy and density. In addition, we report the discovery 
of a putative host galaxy. A comparison of our X-ray and optical 
data may require an appreciable amount of extinction and the 
highest intrinsic neutral hydrogen column density for a short 
GRB to date. Our results have implications for the opening an- 
gle distribution and therefore the observed short GRB rate and 
true energy release. 

Unless otherwise noted, all magnitudes in this paper are 
in the AB system and are corrected for Galactic extinction 
in the direction of the burst using E(B — V) = 0.432 mag 
(Schlegel et al. 1998; Schlafly & Finkbeiner 201 1). We employ 
a standard ACDM cosmology with L1 M = 0.27, Ll\ = 0.73, 
and Hq — 1\ km s” 1 Mpc' 1 . 

2. OBSERVATIONS OF GRB 111020A 
2.1. Swift Observations 

GRB 11 1020 A was detected on 2011 October 20 at 
06:33:49.0 UT by the Burst Alert Telescope (BAT) on board 


5 GRB 090510 also exhibits a post-jet-break-like decay in its optical light 
curve (Nicuesa Guelbenzu et al. 2012). 

6 Please note that Berger (2007) erroneously refers to GRB 061006. 

7 A jet break was reported in the light curve of GRB 090426A (Nicuesa 

Guelbenzu et al. 201 1 1; however, the characteristics of its environment and 
prompt emission are more similar to those of long GRBs (Levesque et al. 
2010; Xin et al. 2011). 
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the Swift satellite (Gehrels et al. 2004; Sakamoto et al. 2011c). 
BAT located the burst at a ground-calculated position of R.A. = 
19 h 08 m 06*9 and deck = -38°01'50'. , 3 (J2000) with 2(1 accu- 
racy (90% containment; Sakamoto et al. 2011b). The X-Ray 
Telescope (XRT) commenced observations of the location of 
the burst at St = 72.8 s (where St is the time after the BAT 
trigger) and detected a fading X-ray source (Section 2.2). The 
UV-Optical Telescope (UVOT) began observations of the field 
at St = 79 s, but no corresponding UV or optical source was 
found within the XRT position. The 3<r limit in the white filter, 
which transmits over X = 1600-8000 A (Poole et al. 2008), 
is 20.3 mag (not corrected for Galactic extinction; Oates & 
Sakamoto 201 1). 

The gamma-ray emission consists of a single pulse with a 
duration of 7go = 0.40 ± 0.09 s in the 15—350 keV band, 
classifying GRB 111020A as a short burst (Sakamoto et al. 
2011b). The spectrum is best fit with a single power law 
with index 1.37 ± 0.26 and a fluence of f Y = (6.5 ± 1.0) x 
I 0“ 8 erg cm -2 (15-150 keV). Spectral lag analysis is not 
conclusive, and there is no clear evidence for extended emission 
(Sakamoto et al. 201 la). 

2.2. X-Ray Observations 

The XRT promptly located a fading, uncataloged X-ray 
source (Evans et al. 2007, 2009; Sakamoto et al. 201 lc) with a 
UVOT-enhanced position of R.A. = 19 h 08 m 12 f !53 and deck = 
— 38°00 , 43'.'8 (J2000) and an uncertainty of 1"6 (Osborne et al. 
2011). XRT observations of the field continued until the source 
faded below the detection threshold at St ~ 3.5 days. 

We also observed the field of GRB 11 1020 A with the 
European Photon Imaging Camera (EPIC-PN) on board the 
X-ray Multi-Mirror Mission ( XMM-Newton ) starting at St = 
0.65 days. With 13.5 ks of on-source observations, we clearly 
detect a source in the energy range 0.5-10 keV, consistent 
with the Swift-X RT position. In addition, we obtained two 
sets of 20 ks observations with the Advanced CCD Imaging 
Spectrometer (ACIS-S, 0.3-10 keV) on board the Chandra 
X-Ray Observatory at St — 2.9 and 10.1 days to refine the 
astrometry and monitor the light-curve evolution. We detect the 
X-ray afterglow in the first Chandra observation but do not 
detect any source at the same location in the second epoch. 

2.2.1. Data Analysis and Spectral Fitting 

We analyze the XRT data using the latest version of the 
HEASOFT package (ver. 6.11) and relevant calibration files. 
We apply standard filtering and screening criteria and generate a 
count rate light curve following the prescriptions from Margutti 
et al. (2010, 2012). Our re-binning scheme ensures a minimum 
signal-to-noise ratio of S/N = 4 for each temporal bin. We 
analyze the XMM data using standard routines in the Scientific 
Analysis System (SAS) ver. 1 1 . We remove the first 5 ks of data 
due to high background contamination, giving a total exposure 
time of 13.5 ks. We extract count rates from a 20" radius 
aperture, and the background is calculated using 20" radius 
source-free regions on the same chip. We use the CIA0 data 
reduction package for the Chandra data. For the first epoch, 
we use a 2" 5 radius source aperture centered on the Chandra 
position and a background annulus with inner and outer radii of 
10" and 35", respectively, giving a source detection significance 
of ~30er. For the second epoch, we extract 1 count in a 2" 5 
aperture at the location of the source, consistent with the average 
background level calculated from source-free regions on the 


2 


The Astrophysical Journal, 756:189 (12pp), 2012 September 10 


Fong et al. 


Table 1 

GRB 1 1 1 020A X-Ray Spectral Fit Parameters 


Telescope 

Detector 

8t 

N ma C b 

pa,b 

C-stat u /dof 



(ks) 

(10 22 cm^ 2 ) 



Swift 

XRT 

0.08-60 

1.0 ±0.3 

2.2 ±0.5 

0.86/188 

XMM 

EPIC-PN 

61.4-76.8 

0 65 +0 - 21 

u. 00—023 

2.0 ±0.4 

1.0/256 

Chandra 

ACIS-S 

250.5-268.5 

0 4 +2 - 3 
u -^-0.4 

1 1 +2J 
1 ‘ — 0.8 

0.32/661 

Swifl+XMM 

XRT+EPIC-PN 

see above 

0 7 <0.20 
u - /o -0.18 

2.0 ±0.3 

0.94/446 


Notes. 

a These values assume a Galactic column density of N\ [ g a ] = 6.9 x I O 20 cm 2 (Kalberla et al. 2005), using an XSPEC model of 
tbabs x ztbabs x pow at z = 0. 
b Uncertainties correspond to a 90% confidence level. 


Table 2 

GRB 1 1 1020 A X-Ray Observations 


8t 

Time Bin Duration 

Unabs. Flux (0.3-10 keV) 

(s) 

(s) 

(erg cm -2 s -1 ) 

6.18 x 10 la 

Swift /XRT 
7.44 x 10° 

(2.38 ±0.79) x 10~ 10 

1.35 x 10 2 

4.64 x 10 1 

(3.80 ± 1.03) x 10“ u 

2.66 x 10 2 

1.71 x 10 2 

(1.64 ±0.42) x 10“ u 

4.15 x 10 2 

1.26 x 10 2 

(2.43 ±0.64) x 10“ u 

5.96 x 10 2 

2.36 x 10 2 

(1.03 ±0.26) x 10“ u 

7.97 x 10 2 

1.66 x 10 2 

(1.83 ±0.49) x 10“ u 

1.14 x 10 3 

5.20 x 10 2 

(8.85 ± 1.72) x 10~ 12 

5.94 x 10 3 

2.46 x 10 3 

(1.65 ±0.33) x 10~ 12 

1.17 x 10 4 

2.46 x 10 3 

(1.07 ±0.27) x 10- 12 

1.94 x 10 4 

6.38 x 10 3 

(9.19 ±2.28) x lO^ 13 

2.58 x 10 4 

6.28 x 10 3 

(1.19 ±0.32) x 10~ 12 

3.19 x 10 4 

5.90 x 10 3 

(1.05 ±0.27) x 10~ 12 

4.29 x 10 4 

1.61 x 10 4 

(8.36 ±2.41) x 10~ 13 

1.26 x 10 5 

1.51 x 10 5 

(1.63 ±0.55) x 10^ 13 

3.09 x 10 5 

2.14 x 10 5 

(1.11 ±0.42) x 10~ 13 

6.91 x 10 4 

XMM/ EPIC-PN 
1.35 x 10 4 

(2.66 ±0.19) x 10^ 13 

2.61 x 10 5 

Chandra/ ACIS-S 
1.98 x 10 4 

(5.96 ±0.89) x 10~ 14 

8.84 x 10 5 

1.98 x 10 4 

<8.95 x 10- 15 


Notes. Upper limits are 3 a . 

a These points were excluded from the broken power-law fit. 


same chip. We therefore take the 3cr background count rate as 
an upper limit on the X-ray afterglow. 

To extract a spectrum from the X-ray data, we fit each of 
the data sets with an absorbed power-law model ( tbabs x 
Ztbabs x pow within the XSPEC routine) characterized by 
a photon index, T, and intrinsic neutral hydrogen absorption 
column, Afn,int> in excess of the Galactic column density in 
the direction of the burst, /Vh.mw = 6.9 x 10 20 cm -2 (typical 
uncertainty of ~ 1 0%; Kalberla et al. 2005; Wakker et al. 
2011), using Cash statistics. For XRT, we utilize data in 
the time interval 0.08-60 ks where there is no evidence for 
spectral evolution. We find an average best-fitting (C-stat„ = 
0.86 for 188 dof) spectrum characterized by T = 2.2 ± 0.5 
and N hm = (1.0 ± 0.3) x 10 22 cm- 2 at z = 0 (Table 1). 
Uncertainties correspond to the 90% confidence level (CL). 
Our best-fit parameters are fully consistent with the automatic 
spectrum fit produced by Mangano & Sakamoto (2011). The 
XMM data are best modeled with a power law characterized 
by T = 2.0 ± 0.4 and yV Rint = (0.65 ± 0.22) x 10 22 cm' 2 
(C-stat„ = 1.0 for 256 dof), consistent with the XRT model 



Figure 1. Unabsorbed X-ray flux light curve for GRB 11 1020 A from Swift- 
XRT (red), XMM (green), and Chandra (blue). Flux errors are la. The Chandra 
3a upper limit is denoted by the blue triangle. The best-fit broken power-law 
model (black solid line) for GRB 111020A is characterized by a\ = —0.78, 
«2 = — 2.1, and tj = 2.0 days. A single power-law model with a = —0.78 
(black dotted) violates the Chandra upper limit. Also plotted are X-ray light 
curves for short GRBs 05 1221 A (dark gray circles; Soderberg et al. 2006; 
Burrows et al. 2006) and 050724 (light gray circles; Grupe et al. 2006). The 
data for GRBs 05 1221 A and 050724 have been scaled for clarity. Gray lines 
trace the afterglow evolution with a break for GRB 05 122 1 A at ~5 days and no 
break for GRB 050724A to ^22 days. 

(A color version of this figure is available in the online journal.) 


parameters. We also fit the first epoch of Chandra data, and 
the resulting parameters are consistent with those from XRT 
and XMM, albeit with large error bars due to low count statistics 
(Table 1). Since we find no evidence for spectral evolution in the 
XRT data, we perform a joint XRT +XMM spectral analysis to 
obtain the best constraints on T and ;V H , m | . The resulting best-fit 
model has T = 2.0 ± 0.3 and Nn,mt = (0.8 ± 0.2) x 10 22 cm~ 2 
(90% CL, C'-stat,, = 0.94 for 446 dof). Although the redshift of 
the burst is unknown, we note that T remains unchanged within 
its la value for z < 3, and we find evidence for intrinsic /V H ,i n i 
in excess of the Galactic value at 6.5a confidence. The best- 
fitting spectral parameters for each of the three data sets and the 
joint fit are summarized in Table 1. 

Applying these parameters to the individual XRT, XMM , 
and Chandra data sets, we calculate the count-rate-to-flux 
conversion factors and hence their absorbed and unabsorbed 
fluxes (Table 2). Applying these conversion factors results in 
the X-ray light curve shown in Figure 1 . 
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hnini-S/GMOS* 
land 20x1 80s 
loch 1 + Epoch 2 


Image subtraction 
Epoch 1 - Epoch 2 


subtraction 


Figure 2. Optical r'-band observations obtained with CMOS on Gemini-South. Left: combined stack of two nights of GMOS i-band data. Stars S 1 and S2 are labeled, 
as well as galaxies G1 and G2. X-ray positions of GRB 1 1 1020A are denoted by the circles (red: Swift-X RT, 1/6 radius, 90% containment; blue: Chandra, 0/33 radius, 
90% confidence). Center: PSF-subtracted image with the centroids of SI and S2 (magenta circles). The subtraction reveals a third source, G3, with i ~ 24.3 mag. 
Right: digital image subtraction of the two epochs obtained at 17.7 hr and 1.7 days, respectively, reveals no residuals in or around the Chandra position. 

(A color version of this figure is available in the online journal.) 


2.2.2. Differential Astrometry 

In the absence of the detection of an optical afterglow 
(Section 2.3), we use our Chandra observations to refine the 
Swift /X RT position (1/6 uncertainty) to subarcsecond accu- 
racy. We perform differential astrometry between our Chandra 
and Gemini Multi-Object Spectrograph (GMOS) observations 
(Section 2.3). To achieve the maximum signal-to-noise ratio, 
we combine both epochs of GMOS i-band observations and 
use SExtractor to determine the positions and centroid un- 
certainty of sources in the field. Performing an absolute as- 
trometric tie to Two Micron All Sky Survey (2MASS) using 
~70 common point sources, we find a resulting rms value of 
<tgmos- 2 mass = O'.' 17 (O'.' 12 in each coordinate). 

To refine the native Chandra astrometry and determine the 
location of the X-ray afterglow relative to the GMOS image, 
we perform differential astrometry. We use CIA0 routines 
mergeall to combine the two Chandra epochs and wavdetect 
to obtain positions and I a centroid uncertainties of X-ray 
sources in the field. We also use wavdetect to determine the 
Chandra position of the afterglow from the first epoch and 
find a la centroid statistical uncertainty a x , ag = O'/ 08. We 
calculate an astrometric tie based on two X-ray and optically 
bright common sources and find weighted mean offsets of 
<5R.A. = -0/27 ± 0/06 and <5decl. = +0/05 ± 0/05, giving 
a tie uncertainty of acxo-GMOS = 0/08. There are three 
additional common, but fainter, sources. An astrometric tie 
using all five sources gives weighted mean offsets and errors 
of SR. A. = —0/29 ± 0/15 and <5decl. = +0/06 ± 0/16, fully 
consistent with our results from using the two bright sources 
alone. We therefore use the astrometric solution from the two 
bright sources only. Applying this solution, we obtain a Chandra 
X-ray afterglow position of R.A. = 1 9 h 08 rn 1 2i49 and deck = 
— 38°00'42/9 (denoted by the blue circle in Figure 2) with a 
total la uncertainty of 0/20, accounting for the 2MASS-GMOS 
astrometric tie, GMOS -Chandra tie, and afterglow positional 
uncertainty. This position is consistent with the XRT position 
and is offset by 0/94 from the XRT centroid. 

2.3. Optical Observations and Putative Host Galaxies 

We initiated i-band observations of GRB 111020A with the 
GMOS mounted on the Gemini-South 8 m telescope on 2011 
October 21.01 UT (St = 17 .7 hr). The data were reduced using 


http ://sextractor. sourceforge . net/ 


the gemini/gmos package in IRAF. In a stack of 9 x 180 s 
exposures in 0/74 seeing and photometric conditions (Figure 2), 
we do not detect any sources within the enhanced XRT error 
circle or coincident with the Chandra position. However, the 
outskirts of the XRT position are partially contaminated by 
emission from a nearby i = 1 7.7 mag star (SI) and a fainter star 
(S2)with/ = 22.7 mag (corrected for A ; = 0.73 mag. Figure 2). 
We detect two additional nearby sources: a faint galaxy (Gl) 
located 2/ 8 away from the center of the Chandra position and a 
brighter galaxy (G2) with a 6/5 offset (Figure 2). 

To search for a fading optical afterglow, we obtained a second, 
deeper set of i-band observations (11 x 180 s) with GMOS at 
St = 1.74 days in 0/67 seeing. Digital image subtraction using 
the ISIS software package (Alard 2000) reveals no variation 
between the two epochs inside or near the X-ray afterglow 
error circles (Figure 2). To calculate the upper limit on the 
afterglow, we add several point sources of varying magnitudes 
in the range i = 24-26 mag around the position in the first 
epoch using IRAF routine addstar. We perform photometry in 
the residual image in 2" apertures using the standard published 
/-band zero point for GMOS-S and place a 3a limit on the 
afterglow of i > 24.4 mag (F v < 0.63 /xJy). We also perform 
photometry in a 1/8 aperture for Gl and a 2/3 aperture for G2, 
giving magnitudes of i = 23.9 ± 0.2 mag and i = 21.9 ± 0.1, 
respectively (Table 3). 

In addition, we obtained /--band observations (3 x 360 s in 
0/62 seeing) with the Low-Resolution Survey Spectrograph 
3 (LDSS3) mounted on the Magellan/Clay 6.5 m Telescope 
concurrent to the first epoch of GMOS observations (St = 
17 .7 hr). We obtained a second, deeper set of observations 
(16 x 150 s in 0/66 seeing) at St fs 180 days, and digital image 
subtraction reveals no evidence for a fading source within the 
X-ray positions to a 3a limit of r > 24.1 mag, where the zero 
point is determined from several standard stars at similar air 
mass. We easily detect G2, with r = 21.9 ± 0.1 mag, but do 
not detect Gl to a 3a limit of r > 24.8 (corrected for A r = 
0.99 mag, Table 3). Finally, we obtained z-band observations 
with LDSS3 (15 x 180 s in 0/60 seeing) at St 180 days. We 
detect Gl at 2.5a significance, z = 24.1 ± 0.4, and G2 with 
z — 21.7 ± 0.1 (corrected for A z = 0.55 mag). 

Since SI and S2 contaminate the Chandra position, we 
subtract their contribution using point-spread function (PSF) 
subtraction on the individual observations and a combined stack 
of the two GMOS epochs. We use standard PSF-fitting routines 
in the IRAF daophot package. Modeling the PSF using four 
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Table 3 

GRB 1 1 1020A Optical Photometry 


Date 

(UT) 

8t 

(days) 

Telescope 

Instrument 

Filter 

Exposures 

(s) 

#FWHM 

(arcsec) 

Afterglow 3 ’ 15 
(AB mag) 

p a,b 

Jjy) 

Gl a 

(AB mag) 

G2 a 

(AB mag) 

G3 a 

(AB mag) 

Aa 

(AB mag) 

2011 Oct 21.01 

0.74 

Magellan/Clay 

LDSS3 

r 

3 x 360 

0.62 

>24.1 

<0.83 

>24.1 

22.00 ± 0.08 

>24.1 

0.987 

2011 Oct 21.01 

0.74 

Gemini- S 

GMOS 

i 

9 x 180 

0.74 

>24.4 

<0.63 




0.734 

2011 Oct 22.01 

1.74 

Gemini-S 

GMOS 

i 

11 x 180 

0.67 






0.734 

2011 Oct 21+22 


Gemini-S 

GMOS 

i 

20 x 180 

0.72 



23.89 ±0.17 

21.91 ±0.05 

24.27 ±0.16 

0.734 

2012 May 17.25 

179.0 

Magellan/Clay 

LDSS3 

z 

15 x 180 

0.60 



24.05 ± 0.41 

21.67 ±0.08 

>23.6 

0.546 

2012 May 17.30 

179.1 

Magellan/Clay 

LDSS3 

r 

16 x 180 

0.66 



>24.8 

21.84 ±0.05 

>24.8 

0.987 


Notes. 

a These values have been corrected for Galactic extinction A, (Schlafly & Finkbeiner 2011). 
b Limits are 3cr . 


bright, unsaturated stars in the field out to a radius of 3" 
(~46Whm) from the center of each star, we subtract several 
stars in the held including SI and S2. The clean subtraction 
of these stars indicates a model PSF representative of the PSF 
of the held. We uncover a faint, mildly extended source (G3) 
on the outskirts of SI at coordinates R.A. = 1 9 h 08 nl 1 2)43 and 
deck = — 38°00'43"07 (J2000). This source, which lies 0"80 
from the center of the Chandra error circle, has a magnitude 
of i = 24.3 ± 0.2 and is a potential host of GRB 111020A 
(Section 3.1). G3 is not detected in the r or z hlters to 3er limits 
of r > 24.8 and z > 23.6 (Table 3). Based on the limited color 
information, r — i > 0.5 and i — z > —0.67, we cannot rule out 
the possibility that this source is a faint star. 

2.4. Radio Observations and Possible Afterglow 

We observed the position of GRB 111020A with the Ex- 
panded Very Large Array 9 (EVLA) beginning on 201 1 October 
20.95 UT (St = 16.1 hr. Program IOC-145) at a mean frequency 
of 5.8 GHz with a total on-source integration time of 65 min- 
utes. We observed 3C 295 and J1937— 1958 for bandpass/flux 
and gain calibration, respectively, and used standard procedures 
in the Astronomical Image Processing System (AIPS; Greisen 
2003) for data calibration and analysis. With the new wideband 
capabilities of the EVLA (Perley et al. 2011), our data have an 
effective bandwidth of ~ 1.5 GHz after excising edge channels 
and data affected by radio frequency interference. The low dec- 
lination of GRB 111020A and the compact D configuration of 
the array at the time of observation caused significant shadow- 
ing and required the removal of 7 out of 27 antennas (the north 
arm of the EVLA). 

Taking into account the highly elongated beam (33" x 7" 
with a position angle of 170°), we detect a 3.7a source with 
an integrated flux density of 48 ± 13 /xJy located at R.A. = 
19 h 08 m 12H0, deck = -38°00'41"2 (5R.A. = l'.'l. Meek = 
3"6, la uncertainty), consistent with the Chandra position. 
The position, peak flux, and integrated flux of the source are 
consistent regardless of our choice of weighting, or if we analyze 
the upper and lower sidebands separately. However, despite the 
statistical significance of the detection, we cannot completely 
rule out sidelobe contribution from nearby bright sources in 
the field due to the low declination of the burst. Therefore, we 
conservatively adopt a 3a upper limit of 39 /jJy on the radio 
afterglow of GRB 1 1 1020A for our analysis. We note that if the 
source is indeed real, then upper limits inferred from the radio 
data can be treated as actual values. 


9 Newly renamed the Karl G. Jansky Very Large Array. 



Figure 3. Probability of chance coincidence, P(<8R), as a function of angular 
distance from the center of the Chandra afterglow position for the three host 
galaxy candidates of GRB 1 1 1020A. The galaxy G3 has the lowest probability 
of chance coincidence, P(<8R ) = 0.01, and is therefore the most probable host 
of GRB 111020A. 


3. RESULTS 

3.1. Galaxy Probabilities of Chance Coincidence 

To assess which galaxy is the most probable host of 
GRB 11 1020 A, we adopt the methodology of Bloom et al. 
(2002) and Berger (2010) to calculate the probability of chance 
coincidence P(< SR) at a given angular separation SR. We de- 
termine the expected number density of galaxies brighter than 
a measured apparent magnitude, m, using the results of deep 
optical galaxy surveys (Hogg et al. 1997; Beckwith et al. 2006): 

a« m) = x io°-33("'-2 4 )-2. 44 arcsec -2 . (1) 

0.33 x ln(10) 

Then the probability of chance coincidence is given by (Bloom 
et al. 2002) 

P(< SR) = 1 - (2) 

We calculate P(<8R) for each of the three host galaxy 
candidates (Figure 3) and find that G3 is the most probable 
host of GRB 1 1 1020 A with P(<SR) = 0.01, while for G1 and 
G2 the values are P(<SR ) = 0.10 and 0.12, respectively. 
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3.2. X-Ray Light-curve Fitting and a Jet Break 


The temporal behavior of the X-ray afterglow flux is charac- 
terized by a steady power-law decline until St ^ 2 days, when 
there is a significant steepening in the light curve (Figure 1). A 
single power-law model with a decline rate determined by the 
X-ray data at early times (t < 2 days) provides a poor fit to the 
late-time data (dotted line in Figure 1); in particular, it overesti- 
mates the Chandra detection and upper limit. To quantitatively 
assess the shape of the X-ray light curve, we therefore invoke a 
broken power-law model, given by 


Fx = F x , o 



i A 


(3) 


where F x . o = 2 V'F x (t = tj), a x . i and a x ,2 are the power- 
law indices pre- and post-break, respectively, tj is the break 
time in seconds, and a is a dimensionless smoothness parameter 
that characterizes the sharpness of the break. We perform a 
three-parameter / 2 -grid search over Fx, o, a x. i and tj. If we 
use a relatively sharp break (e.g., s ~ —10), the Chandra 
3a upper limit constrains a x ,2 < — 1.7. If we allow for a 
smoother break (e.g., s ~ —1), a x . i remains unchanged but 
the break occurs at later times ( tj *=» 4 days) and ax , 2 is 
required to have a steeper value of <—2.2 to accommodate 
the Chandra upper limit. This scenario generally provides 
a poorer fit to the last Chandra and Swift/XKT points. We 
therefore adopt the sharp-break scenario. Fixing s = —10 
and ax , 2 = —2.1, we find a best-fit broken power-law model 
characterized by F x (tj) = (1.36 ± 0.45) x 10 -13 erg cm~ 2 s _1 , 
ax,i = —0.78 ± 0.05, and tj = 2.0 ± 0.5 days (la, xl= 1.1 
with 15 dof, Figure 1). This best-fit model is shown in Figure 1. 
The best-fit parameters are independent of our choice of ax, 2 
between —1.7 and —3. We also note the presence of a slight flux 
enhancement relative to the power-law decay at ~ 3 x 1 0 4 s 
(0.35 days). If we remove these points from our fits, the resulting 
best-fit parameters remain unaffected. 

The required change in the temporal index is Aa^ > 0.9. 
There are several possibilities that can explain breaks in GRB 
afterglow light curves. The first scenario is the transition of 
the cooling frequency across the band, but this only predicts 
Aa = 0.25 (Sari et al. 1998). An alternative possibility is the 
cessation of energy injection, from either refreshed shocks or a 
long-lasting central engine (e.g., Rees & Meszaros 1998; Sari 
& Meszaros 2000; Zhang & Meszaros 2002), which has been 
used to explain the termination of a shallow decay or plateau 
phase in the X-ray and optical light curves of several long 
GRBs. However, these cases all exhibit earlier temporal breaks 
at ~10 3 -10 4 s with Aa 12 ~ 0.7 (ax,i ~ —0.5, ax , 2 ~ —1-2; 
Nousek et al. 2006; Zhang et al. 2006; Liang et al. 2007). 
Attributing the break in GRB 111020A to the cessation of 
central engine activity would require sustained energy injection 
from the start of XRT monitoring to the break time, ~100 s 
to 2 days, whereas the timescales of energy injection for long 
GRBs are < few hours (Nousek et al. 2006; Zhang et al. 2006; 
Liang et al. 2007; Racusin et al. 2009). Single episodes of 
energy injection have also been observed in two short GRBs: 
051221A and 050724A (Berger et al. 2005; Soderberg et al. 
2006; Burrows et al. 2006; Grape et al. 2006)). The light 
curve of GRB 05 1221 A, which exhibits a power-law decay 
with index a x l = —1.1, a plateau, and a return to the same 
power law (Aai 2 = 0), is interpreted as a single period of 
energy injection (Soderberg et al. 2006; Burrows et al. 2006). A 
superimposed flare on the light curve of GRB 050724A with a 


single underlying decay index of a x 1 = —0.98 is also possibly 
related to late-time reactivation of the central engine (Berger 
et al. 2005; Grape et al. 2006; Figure 1). Neither of these light 
curves resembles the behavior of GRB 111020A, where the 
change in slope is substantially greater. 

Yet another possibility to explain the break is a sharp change 
in the external density. However, models for density jumps in a 
uniform medium (Nakar & Granot 2007) predict that the density 
would need to decrease by greater than a factor of ~10 3 to 
account for the observed Aa^ > 0.9 steepening. More realistic 
density contrasts of ~10 predict Aa max 0.4 in optical and 
X-ray afterglow light curves (Nakar & Granot 2007). 

Finally, we consider that the observed steepening is a jet 
break, when the edge of a relativistic ally beamed outflow 
becomes visible to the observer and the jet spreads laterally 
(Sari et al. 1999; Rhoads 1999). This model is often adopted to 
explain Aa^ ~ 1 in the light curves of long GRBs (e.g., Frail 
et al. 2001 ; Bloom et al. 2003; Racusin et al. 2009) and has been 
observed in one other short burst, GRB 051221A (A an ~ 0.9, 
Figure 1; Soderberg et al. 2006; Burrows et al. 2006). Given the 
similarity in Aan and the timescales of jet breaks in both short 
and long GRBs, we conclude that the observed steepening in 
the light curve of GRB 1 1 1020A is best explained by a jet break 
at tj = 2.0 ± 0.5 days. 

3.3. Afterglow Properties 

We utilize our radio, optical, and X-ray observations to con- 
strain the explosion properties and circumburst environment of 
GRB 1 1 1020A. In particular, we adopt the standard synchrotron 
model for GRB afterglows (Sari et al. 1998; Granot & Sari 
2002), which provides a mapping from observable properties 
to the isotropic-equivalent kinetic energy (.Ejc.iso), circumburst 
density (no), and the fractions of post-shock energy in radiating 
electrons ( e e ) and magnetic fields (e«). We use data at the time 
of the radio and first optical observations (St = 17.7 hr), as well 
as the decay indices from the full X-ray light curve. 

First, we constrain the electron power-law index p, using a 
combination of temporal and spectral information. From the 
X-ray light curve, we measure ax , 2 < —1.7 (Section 3.2). For 
p = —ax, 2 , appropriate for a spreading jet (Sari et al. 1999), 
we can then constrain p > 1.7. To further constrain p and 
investigate the location of the cooling frequency, v c , we com- 
pare the values ax,i = —0.78 ± 0.05 and j3 x = —1.04 ± 
0.16 (f x = I — T ; la) to the closure relations for a rela- 
tivistic blast wave in a constant-density interstellar-medium- 
like medium for p > 2, a typical environment expected for a 
short GRB from a non-massive star progenitor (Sari et al. 1999; 
Granot & Sari 2002). If v c > v x , then the independently derived 
values for p from the temporal and spectral indices are incon- 
sistent: p = 2.0 ± 0.07 from ax.i, and p — 3.1 ±0.32 from f x 
(errors are la). 

However, if v c < v x , we obtain p — 1.7 ± 0.07 from ax.i 
(Granot & Sari 2002), which is consistent with the p value 
inferred from a X 2 but yields a divergent total integrated energy 
in electrons unless a break at high energies in the distribution is 
invoked. Although a flat electron distribution (p < 2) is possible 
and not uncommon (e.g., Dai & Cheng 2001; Panaitescu & 
Kumar 2001; Racusin et al. 2009), the standard relations for 
1 < p < 2 yield p = 0.84 ± 0.25 from ax.i - This solution is 
not self-consistent and would also require an unusually sharp 
break of Ap > 1.2 in the electron distribution. Continuing 
with the assumptions that v c < v x and p > 2, we obtain 
p — 2.1 ± 0.32 from /3 X , which is marginally consistent with 
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Itf* 10 -3 10' 2 3 10 _1 10° 10“ ICf 3 10" 2 3 10" 1 10° 

n 0 (cm- 3 ) n Q (cm" 3 ) 

Figure 4. Isotropic-equivalent kinetic energy and circumburst density parameter space for GRB 11 1020 A assuming € e = €b = 1/3 (left) and e e = e# = 0.1 
(right). The lower limit on the density (dark blue) is set by the condition that v c < vx (Equation (6)), while the upper limit (light blue) is set by radio observations 
(Equation (5)). Also plotted are the values for 2?K,iso at z = 0.1, 1, and 3 (red). The allowable parameter space set by these constraints is tilled in yellow. 

(A color version of this figure is available in the online journal.) 


the value inferred from the temporal index. Put another way, 
a — 3/3/2 — 0.77 ± 0.30, which satisfies the closure relation for 
v c < Vx (Sari et al. 1998). We therefore conclude that v c < v x . 
We note that the spectral index is generally more reliable in the 
determination of p because it remains constant over time and 
is not subject to alternative processes such as energy injection 
or flaring. In this case, the same fi x was also independently 
determined from both the XMM and XRT data sets (Table 1). 
Therefore, for the rest of our calculations, we take a reasonable 
value of p = 2.1 as determined from f J > x . 

We next determine a set of constraints on no and T K ,iso based 
on the X-ray flux density, radio limit, and the condition that 
v c < v x . First, we use the X-ray afterglow emission as a proxy 
for TiK.iso assuming that the X-ray emission is from the forward 
shock. For v c < v x at the time of our broadband observations 
( St = 17.7 hr), we use F x = 0.032 jxJy at v x = 2.4 x 10 17 Hz 
(1 keV) and p = 2.1 to obtain (Granot & Sari 2002) 

E KMo « 2.2 x 10 50 (1 + zrV'°V' 02 X ; 28 (4) 

where d\_p% is the luminosity distance in units of 10 28 cm. 
Next, we use £k i S0 to constrain «o- Using our 3er EVLA 
limit of /Vrad < 39 //Jy, we can determine an upper limit on 
no under the reasonable assumption that our observed radio 
band (v = 5.8 GHz) is above the self-absorption frequency 
(v a < v ra d < v m , Fj, ra d oc n Q ) at the time of observations. For 
this scenario (Granot & Sari 2002), 

n 0 < 1.7 x 10" 3 £" 5 / o 3 52 (1 +z)- 5/3 ^ 28 e e 4/3 e" 2/3 cm” 3 , (5) 

where T K ,iso.52 is in units of 10 52 erg and no is in cm” 3 . As noted 
in Section 2.4, if the marginal radio detection is indeed real, then 
this upper bound can be replaced with an equality. Finally, we 
can place a lower limit on the density using the condition that 
v c < v x (i.e., v c < 2.4 x 10 16 Hz; 0.1 keV): 

n 0 > 4.5 x 10” 4 (1 + z)~ l,2 e B 3/2 E~^l 52 cm” 3 . (6) 

It is clear that T Kis „ and n () depend sensitively on our choice 
of z,€ e , and €b- The fractions e e , eg are not expected to exceed 
~l/3. We therefore calculate E K , so for two representative 
cases: e e = e B = 1/3 (Case I), and more typical values of 


e e = e B = 0. 1 (Case II). We then calculate the range of allowed 
no set by Equations (5) and (6), which becomes wider as the 
redshift increases. 10 For Case I, this requires that z. > 0.2, 
below which the constraints on n q conflict (Figure 4). At the 
median observed redshift of the short GRB population, z ~ 0.5, 
we obtain E^ xso ~ 3.7 x 10 50 erg and no = 0.01-0.06 cm” 3 . 
For Case II, the constraints on no require a higher redshift of 
Z > 1.5 (Figure 4). For a fiducial redshift of z = 1.5, we obtain 
^K.iso ^ 1.2 x 10 52 erg and no = 0.008 cm” 3 . The parameters 
for the two cases are listed in Table 4. Although we cannot 
distinguish between these two scenarios, both cases require low 
circumburst densities of n ~ 0.01-0.1 cm” 3 . 

3.4. Jet Opening Angle 

In the context of a jet break, we use the time of the break from 
the X-ray light curve (2.0±0.5 days) and the circumburst density 
and energy estimates from the previous section to constrain 0j . 
The time of the break is a direct reflection of the jet opening 
angle (Sari et al. 1999; Frail et al. 2001), 

6j = 0. lr 3/ d 8 (l + z)- 3/8 £k,/o, 52 «/ 8 ’ (7) 

where tj o is expressed in days. For our fiducial Case I (z = 
0.5, € e = €b = 1/3), fiK.iso ~ 3.7 x 10 50 erg and n rs 
0.01-0.06 cm” 3 give 0j = 7°-8°. This leads to a beaming 
correction on the energy of /* = [1 — cos(0,-)] = 0.007-0.01 
and therefore a true kinetic energy E x = //£K,iso ^ (3-4) x 
10 48 erg (Table 4). To estimate the beaming-corrected y- ray 
energy, we infer E y ; so from the Swift/ BAT fluence and apply 
a bolometric correction factor of five to roughly convert to a 
representative observed y-ray energy range of ~ 10-1000 keV. 
This factor is derived from short GRBs observed by satellites 
with wider energy coverage (Berger 2010; Margutti et al. 2012). 
We obtain E Yt ; so = 2.1 x 10 50 erg and therefore a true y-ray 
energy of E r 2 x 10 48 erg. 

For Case II (z = 1.5, e e = £b = 0-1), where «o ^ 
0.008 cm” 3 and Fk.iso ~ I -2 x 10 52 erg, we obtain a smaller 


10 Assuming non-equipartition in two alternative cases, 65 = 0.01, e e =0.1 
and €b = 0.1, e e = 0.01, our constraints on «o conflict unless z 3 and result 
in 720 ~ 0.01-1 cm -3 . Based on our observations of the putative host galaxies 
and the typical luminosity of short GRB hosts, we find a high-z origin unlikely 
(see Section 4.2). 
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Table 4 

Physical Parameters of GRB 1 1 1020A 


Parameter 

Case I [z = 0.5, € e = cb = 1/3] 

Case II [z = 1.5, e e = eg = 0.1] 

tj 

2.0 ± 0.5 days 2 

2.0 ± 0.5 days 2 

Ey, iso 

2.1 x 10 50 erg 

1.9 x 10 51 erg 

£K,iso 

3.7 x 10 50 erg 

1.2 x 10 52 erg 

no 

0.01-0.06 cm' 3 

0.008 cm' 3 

0j 

-~4 

0 

1 

00 

o 

3° 

fb 

0.007-0.01 

0.001 

Ey 

2 x 10 48 erg 

3 x 10 48 erg 

E k 

(3-4) x 10 48 erg 

2 x 10 49 erg 

Etot 

(5-6) x 10 48 erg 

2 x 10 49 erg 

Vy 

0.3-0.4 

0.15 


Note. a Uncertainties correspond to a lcr confidence level. 


opening angle of 0j 3°. This leads to //, fs 1.4 x 10' 3 and 
hence E y 3 x 10 48 erg and E K ~ 2 x 10 49 erg. 

In both cases, the true y-ray energy is a few x 10 48 erg, while 
the kinetic energy is an order of magnitude higher at z = 1.5 than 
at z = 0.5. This results in a lower y-ray conversion efficiency 
(rjy = Ey/Etot ) for Case II of 0.15 compared to 0.3-0. 4 for 
Case I (Table 4). The total energy even for Case II is ~ 10-100 
times lower that for long GRBs. 

3.5. Extinction 

We investigate the presence of extinction by comparing 
the unabsorbed X-ray flux and the optical non-detection at 
St = 17.7 hr. Since we do not know the exact location of 
the cooling frequency, we assume a maximum value v c , max 
of 2.4 x 10 17 Hz (1 keV) and extrapolate the X-ray flux 
to the optical band using the shallowest possible slope of 
ft = —(p — l)/2 = —0.55 to obtain the lowest bound on 
the expected optical afterglow flux in the absence of extinction; 
any other assumption for the location of v c < v x would result 
in a higher predicted optical flux density. For p = 2.1 we 
estimate F vopt ~ 1.1 /xJy (i = 23.8 mag). Given that our 
observed 3er upper limit is i > 24.4 mag, this implies a lower 
limit on the optical extinction in excess of the Galactic value of 
A i > 0.6 mag. 1 1 In the rest frame of the burst for a Milky Way 
extinction curve, this translates to A y ost > 0.6 mag for z = 0.5 
and Ay OSt > 0.2 at z = 1 .5 (Cardelli et al. 1989). Using Galactic 
relations between Nh and Ay, /Vnint/Ay ~ (1.7-2. 2) x 10 21 
(Predehl & Schmitt 1995; Watson 2011), we infer lower limits 
°f A^n.int 10 21 cnT 2 at z — 0.5 and (V H ,mt 4.4 x 10 2 ° cm -2 
at z = 1.5, consistent with our value of 7.5 x 10 21 cm' 2 
(z = 0) derived from the X-ray spectrum (Table 1). However, an 
appreciable extinction is unexpected given the burst’s location 
on the outskirts of its potential host galaxy. We note that the 
burst is located at Galactic coordinates (/, b) — (359?3, —19? 4), 
which is toward the Galactic bulge on a steep gradient in the 
dust map (Schlegel et al. 1998) and thus may be subject to 
substantial (~30%) uncertainties in the Galactic extinction. 12 
Taking this uncertainty into account reduces the required A y 0St 
to >0.2-0. 3 mag depending on the redshift of the burst. 


1 1 We note that for p < 1 .9, no host galaxy extinction is required. 

12 Using a high-resolution (#fwhm = 15 ff ) WISE 12 // m map, we do not see 
strong evidence for any thin-dust filaments at the location of the burst that 
would result in >30% uncertainties in the Galactic Ay (D. Finkbeiner 2012, 
private communication). 


4. DISCUSSION 
4.1. Environment 

From our broadband observations, we constrain the circum- 
burst density of GRB 111020Aton 0 ~ 0.01-0.1 cm -3 , which is 
consistent with the low values inferred for a few previous short 
GRBs (Soderberg et al. 2006; Panaitescu 2006; Stratta et al. 
2007; Perley et al. 2009b; Berger 2010; Fong et al. 2011). The 
inferred density fits well with the framework of NS-NS/NS-BH 
binary progenitor systems, which may be subject to substantial 
kicks from their host galaxies and are predicted to typically oc- 
cur at densities of ~10' 6 to 1 cm' 3 (Perna & Belczynski 2002; 
Belczynski et al. 2006). 

GRB 111020A has an offset of ~(y.'80 from its most prob- 
able host galaxy (G3, Figure 2). For redshifts between z = 
0.5 and 1.5, this translates to a projected physical offset of 
~5-7 kpc, which is comparable to the median value of ~5 kpc 
for well-localized short GRBs with host associations (Fong et al. 
2010; Church et al. 2011). Although G3 has the lowest prob- 
ability of chance coincidence by an order of magnitude (Fig- 
ure 3), we cannot rule out the possibility that G3 is a faint 
star. The next most probable galaxies, G1 and G2, are situ- 
ated 2'.'8 (17-24 kpc) and 6'.'5 (40-56 kpc), respectively, from 
GRB 11 1020 A (Figure 2). If the burst originated from one of 
these galaxies, this would put GRB 111020A at the high end 
of the observed offset distribution, similar to the growing sub- 
class of apparently “hostless” short GRBs, which likely occur 
>30 kpc from their host galaxies (Berger 2010). All of these in- 
ferred offsets are consistent with predicted offset distributions of 
NS-NS/NS-BH binaries originating in Milky- Way-type galax- 
ies (Bloom et al. 1999; Fryer et al. 1999; Belczynski et al. 2006; 
Salvaterra et al. 2010). 

Most short GRB host galaxies with confirmed spectroscopic 
redshifts have measured luminosities of Lb ^ 0.1-1L* (Berger 
et al. 2007). The apparent magnitude of G3 is i ta 24.3, which 
corresponds to L b fs 0.1-1L* overz rs 0.5-2. 3 when compared 
to the luminosity function of galaxies at corresponding redshifts 
in the DEEP2 and LBG surveys (Willmer et al. 2006; Reddy & 
Steidel 2009). This is consistent with the redshift range inferred 
from the afterglow. 

We next investigate the nature of the dust and gas in the 
environment of GRB 1 1 1020A through an analysis of Ay 0st and 
An.int- We have shown that the burst requires dust extinction 
of Ay OSt > 0.2-0. 6 mag, depending on the redshift of the 
burst and the uncertainty in Galactic extinction. We have also 
measured a neutral hydrogen column density intrinsic to the 
burst environment of IVn.int = (7.5 ± 2.0) x 10 21 cm' 2 at z = 0, 
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Figure 5. Excess neutral hydrogen column density, IVH.int* vs. redshift for 
GRB 111020A (lfr, 2a, and 3cr intervals denoted by green, red, and orange 
lines) along with six Swift short GRBs with measured redshifts and optical 
afterglows (black filled circles) and two (GRBs 060801 and 101219A) with 
only X-ray afterglows (open circles). Also plotted are 1 1 short GRBs without 
secure redshifts (gray lines), 4 of which have only upper limits on (gray 

dashed). For GRBs without redshifts, the ,Vn i nt value at z = 0 is scaled by 
(1 + z) 2 ' 6 (Galama & Wijers 2001). Errors and upper limits are at the 90% 
confidence level. The weighted mean for all short GRBs (black line) over the 
redshift interval z = 0-1 .2 is also shown. GRB 1 1 1020A has the highest (Vn.mt 
of a short GRB to date and is well above the mean for short GRBs. 

(A color version of this figure is available in the online journal.) 

which becomes higher for any other choice of z. High values of 
both dust extinction and X-ray absorption have been linked to 
“dark” GRBs (e.g., Perley et al. 2009a; Campana et al. 2012), 
which have optically sub-luminous afterglows compared to their 
X-ray or near-lR (NIR) counterparts and can quantitatively be 
classified by \Pox\ < \fix \ — 0.5 (van der Horst et al. 2009; see 
also Jakobsson et al. 2004). With \fx\ = 1-0 and \fox\ < 0.46, 
GRB 1 1 1020A is consistent with this definition of dark GRBs. 
While optical extinction intrinsic to long GRB environments is 
not uncommon and commensurate with their origin in dusty, 
star-forming regions, evidence for substantial extinction has 
been reported for only one other short burst, GRB 070724 A, 
which required Ay ost > 2 mag to explain the suppression of 
optical emission relative to the NIR (Berger et al. 2009; Kocevski 
et al. 2010). The location of GRB 070724A on the outskirts 
of its host galaxy, ~5 kpc from the center, suggested either 
an origin in a star-forming region or a progenitor system that 
produced the dust itself (Berger et al. 2009). The potentially 
appreciable extinction and the location with respect to its 
putative host suggest that the same conclusions may be drawn 
for GRB 111020A. 

On the other hand, the relation between /V H ,int. and the 
darkness of a burst is less clear. A recent study of long, dark 
GRBs shows them to have higher intrinsic column densities than 
non-dark GRBs, which suggests that the darkness of a burst is 
largely due to absorption by circumburst material (Campana 
et al. 2012). To investigate this relationship for GRB 1 1 1020A, 
we extract spectra and best-fitting An.int for all short GRBs with 
XRT-detected afterglows in the same manner as GRB 1 1 1020A 
(see Section 2.2.1), over time ranges with no evidence for 
spectral evolution. There are 22 short bursts with sufficient 
X-ray counts to perform spectral analysis, 11 of which have 
known redshifts (Table 5). We find a short GRB weighted 
average of Nu- mt (z = 0) = (1.1 ± 0.14) x 10 21 cm -2 (90% 
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Table 5 

Intrinsic X-Ray Column Density of Hydrogen, A^ino for Swift Short GRBs 


GRB 

z 

Wl.mt 

(10 21 cm" 2 ) 

g above Zero 

050724 

0.258 

3 20 +0 97 

5.7 

051210 


<0.54 


051221 A 

0.547 

1 Q? +0 - 73 

4.5 

060313 


0 45+0-36 

2.1 

060801 

1.131 

3 02 +2 22 

2.4 

061006 

0.438 

<2.04 


061201 


n Q4+O.6O 
— 0.53 

2.7 

0707 14B 

0.923 

3.892 1 1 8 6 7 1 

4.2 

070724 A 

0.457 

<1.89 


071227 

0.383 

2 84 +0 ’ 72 
z -°^-0.65 

6.8 

080123 


i i O+0.28 
i - lz -0.26 

6.8 

080905A 

0.122 

2 04 +1-58 

z - uq --1.33 

2.3 

090510 

0.903 

<0.80 


090515 


0 56 +0 - 30 

u.JO—0 27 

3.2 

090607 


<0.79 


091109B 


<1.58 


1001 17A 

0.915 

4 10 +341 

^f.iu_2 7i 

2.2 

100702 A 


4 07+3.67 
^• J/ -3.05 

2.1 

101219A 

0.718 

6.6C 323 

3.3 

110112A 


<0.92 


11 1020 A 


7-50! 2 1 ° 8 

6.5 

111117A 


1-84-i 2 05 

2.6 

111121A 


2.4L°o 8 7 2 4 

5.1 


Note. Errors and upper limits quoted correspond to a 90% confidence level; 
Z = 0 is assumed when the redshift is not known. 


CL, Figure 5). In comparison, GRB 1 1 1020A has a high value 
of V H ,int = (7.5 ± 2.0) x 10 21 cm -2 at z = 0 (Figure 5). Taken at 
face value, it is surprising to find a large (V H , i n t for a substantial 
offset and may suggest that the burst occurred in a relatively 
metal-rich environment. 

4.2. Beaming, Energetics, and Rates 

We uncover a break in the X-ray light curve of GRB 1 1 1020A 
at «s2 days, which we interpret as a jet break (Section 3.2). 
Depending on our values for z, e e , and e#, we infer an opening 
angle of ^3°-8°. This is reminiscent of the first jet break 
discovery in GRB 05 1221 A, with 9j ~ 7° (Soderberg et al. 
2006; Burrows et al. 2006), and suggests that at least a fraction of 
these events are highly collimated. In addition, temporal breaks 
at tj < few hours in GRBs 061201 (Stratta et al. 2007) and 
090510 (De Pasquale et al. 2010; Nicuesa Guelbenzu et al. 
2012), if interpreted as jet breaks, lead to 9j 1° (Figure 6). 
However, these two cases resemble early breaks in long GRBs 
that are generally attributed to the cessation of energy injection, 
and not collimation. 

Although the remaining short GRB afterglow data are sparse, 
the lack of observed jet breaks in their X-ray light curves can 
be used to place lower limits on the opening angles. Indeed, 
Chandra observations of GRB 050724A out to 22 days indicated 
9j > 25°, consistent with a spherical explosion (Grupe et al. 
2006). A recent study by Coward et al. (2012) analyzed the 
sample of short GRB Swift /XKT light curves up to 201 1 August 
with monitoring >1 day, which included six additional events 
and inferred 9j > 6°-16°, assuming «o = 1 cm -3 for all bursts. 
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Figure 6. Distribution of opening angles for long (orange) and short (blue) 
GRBs. Arrows represent upper and lower limits. The long GRB population 
includes pre-Swift (Frail et al. 2001; Berger et al. 2003; Bloom et al. 2003; 
Ghirlanda et al. 2004; Friedman & Bloom 2005), Swift (Racusin et al. 2009; 
Filgas et al. 2011), and Fermi (Cenko et al. 2010, 201 1; Goldstein et al. 2011) 
bursts. The opening angle for GRB 1 1 1020A ranges from ~3° to 8° (depending 
on the redshift), while GRB 05 1221 A has Qj & 7° (Soderberg et al. 2006; 
Burrows et al. 2006). Tentative jet breaks (blue dashed) for GRBs 061201 
(Stratta et al. 2007) and 090510 (De Pasquale et al. 2010; Nicuesa Guelbenzu 
et al. 2012) are at ~1°. Short GRB lower limits are from the non-detection of 
jet breaks in Swift /X RT data (this work; revised from Coward et al. 2012) and 
Chandra data for GRB 050724A (Berger et al. 2005; Grupe et al. 2006). 

(A color version of this figure is available in the online journal.) 


We revise this analysis for three of the events with robust X-ray 
light curves (GRBs 070714B, 070724A, 071227; data analysis 
prescriptions from Margutti et al. 2012) employing a more 
representative uq ~ 1 0 2 cm -3 (e.g., Soderberg et al. 2006 and 
this work). We derive E y - VM from the reported fluences, applying 
a bolometric correction when necessary to represent an energy 
range of ~1 0-1 000 keV, and infer more realistic lower limits 
of >2°-6° (Figure 6). These limits are indeed lower than the 
detected values for GRB s051221A and 111 020 A and therefore 
do not add strong constraints on the distribution. We caution 
that the sample presented here represents only the ~30% of the 
Swift short GRB population that has bright X-ray afterglows and 
relatively slow flux decline rates; the remaining fraction does 
not have detectable X-ray afterglows or fades too quickly, so 
constraints cannot be placed on their collimation. 

There are now two short GRBs with opening angle measure- 
ments, two with measurements based on more tentative early 
breaks, and an additional four with lower limits (Figure 6). 
These early constraints create a distribution that may mimic the 
distribution for long GRBs, which ranges from ~2° to 20° with 
a median of 7° (Figure 6; Frail et al. 2001; Berger et al. 2003; 
Bloom et al. 2003; Ghirlanda et al. 2004; Friedman & Bloom 
2005; Racusin et al. 2009; Cenko et al. 2010; Filgas et al. 2011; 
Goldstein et al. 2011; Cenko et al. 2011). More events are needed 
to assess the real differences between the distributions. How- 
ever, simulations of post-merger black hole accretion predict 
jets with 0j ~ 5°-20° (Aloy et al. 2005; Rosswog 2005; Rez- 
zolla et al. 2011) to several tens of degrees (Ruffert & Janka 
1999b; Rezzolla et al. 2011) depending on the mechanism of 
energy extraction and Lorentz factor, so there are expectations 
on theoretical grounds that the short GRB distribution is wider. 
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The first major ramification of collimation is the correc- 
tion to the total energy release: the true energy is lower than 
the isotropic-equivalent value by the beaming factor, f,. For 
GRB 111020A with an opening angle of »3°-8°, this cor- 
rection factor is substantial, 0.001-0.01. Depending on the 
redshift, the beaming-corrected energy of GRB 11 1020 A is 
E y (2-3) x 10 48 erg (Table 4), which is an order of magni- 
tude lower than for GRB 05 1221 A, with E y ~ (1-2) x 10 49 erg 
(Soderberg et al. 2006; Burrows et al. 2006), and GRB 050724A, 
with E y R* (0.4-4) x 10 5 ° erg (Grupe et al. 2006). The 
three remaining events with opening angle lower limits, 
GRBs 070714B, 070724A, and 071227, have ranges of E y % 
10 48 -10 51 erg, where the upper bound is set by the isotropic- 
equivalent y-ray energy in the 0-1 000 keV band. The small 

population of short GRBs with measured E y therefore has a 
median value of E y ~ 10 49 erg, which is an order of magnitude 
below Swift long GRBs (Kocevski & Butler 2008; Racusin et al. 
2009) and two orders of magnitude below the pr e-Swift popu- 
lation (Frail et al. 2001; Bloom et al. 2003). Again, this sample 
is incomplete because we can only measure E Y for bursts with 
well-constrained opening angles. 

In a similar vein, we compare the beaming-corrected kinetic 
energy and total energy ( Ek , E lM ) of GRB 111020A to the 
values for other short bursts. Because ZtK.iso is more sensitive 
to our choices for z, € e , and €b, we infer different values for 
Cases I and II. For Case I, we infer Ek R» (3-4) x 10 48 erg, 
£ tot = Ey + E k R» (5-6) x 10 48 erg, and rj Y R» 0. 3-0.4. For 
Case II, we calculate Ek ~ 2 x 10 49 erg, Zs tot » 2 x 10 49 erg, 
and rjy R* 0.15. GRB 051221A had E K ™ 8 x 10 48 erg and a 
total energy release of ~2.5 x 10 49 erg (Soderberg et al. 2006; 
Burrows et al. 2006), while GRB 050724 had a total energy of 
10 50 — 10 51 erg. With E tot % (0.5-2) x 10 49 erg, GRB 111020A 
may be on the low end of the total energy distribution, but more 
events with beaming-corrected energies are needed to better 
characterize the distribution for short GRBs. 

The true total energy release of short GRBs has strong 
implications on the energy extraction mechanism. Two primary 
mechanisms, the thermal energy release from vv annihilation 
in a baryonic outflow (laroszynski 1993; Mochkovitch et al. 
1993) and MHD processes in the black hole’s accretion remnant 
(e.g., Blandford & Znajek 1977; Rosswog et al. 2003), give 
different estimates for the expected energy release. Predictions 
for vv annihilation are largely dependent on the mass of 
the disk and efficiency to produce pairs. Simulations of an 
outflow due to vv annihilation suggest that beaming-corrected 
total energy releases could reach 10 48 — 10 49 erg (Ruffert & 
lanka 1999a, 1999b; Popham et al. 1999; Rosswog 2005; 
Birkl et al. 2007; Lee & Ramirez-Ruiz 2007). Higher energy 
releases can be obtained from MHD processes, which can 
produce luminosities of >10 52 erg s -1 (>10 50 erg s -1 when 
corrected for beaming; Popham et al. 1999; Rosswog et al. 
2003; Lee & Ramirez-Ruiz 2007) depending on the nature of 
the magnetic field amplification. While the true energy releases 
of GRBs 05 1221 A and 050724A suggest that MHD processes 
may be powering these events (Berger et al. 2005; Grupe et al. 
2006; Soderberg et al. 2006; Burrows et al. 2006), the total 
energy of GRB 1 1 1020A is consistent with predictions for both 
scenarios. 

The second major consequence of beaming is that the true 
event rate is higher than the observed rate by the inverse of 
the beaming factor (i.e., R tme = fjf 1 R 0 bs)- Thus, beaming 
provides essential information for understanding the relation 
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to various progenitor systems and is of particular interest since 
the NS-NS/NS-BH merger rates, which are a critical input for 
estimates of Advanced LIGO gravitational wave detections, are 
highly uncertain (e.g., Abadie et al. 2010; Metzger & Berger 
2012). The current estimated observed short GRB volumetric 
rate is ~10 Gpc' yr -1 (Nakar et al. 2006). The estimated 
NS-NS merger rate is much higher: ~200-3000 Gpc -3 yr“ 1 
(Kalogera et al. 2004; Nakar et al. 2006). 

The discrepancy in these rates can be explained if short 
GRBs have typical 0j ~ 8° (fjf 1 ~ 100; see also Metzger 
& Berger 2012). The determination of GRB 111020A’s open- 
ing angle of 3°-8° (ff 1 = 100-730), along with the small 
but increasing sample of opening-angle constraints for short 
GRBs, implies that at least a fraction of these events are sig- 
nificantly beamed and that the true rate of short GRBs is at 
least ~ 1 00-1000 Gpc -3 yr“ 1 . With a few additional opening- 
angle measurements, this value can be significantly improved. 
Other proposed progenitor models, e.g., WD-WD mergers or 
accretion-induced collapse of a WD/NS (Qin et al. 1998; Levan 
et al. 2006; Metzger et al. 2008), have estimated rates of 
<1000 Gpc 3 yr” 1 and ~0. 1-100 Gpc -3 yr~ 1 , respectively (Lee 
& Ramirez-Ruiz 2007; Darbha et al. 2010), so if a large fraction 
of short GRBs have opening angles of <25°, these systems may 
not contribute significantly to the progenitor population. 

5. CONCLUSIONS AND FUTURE WORK 

We have presented observations of GRB 1 1 1020A, utilizing 
extensive coverage in the X-rays with Swift /X. RT, XMM, and 
Chandra to uncover a temporal break, most naturally explained 
as a jet break. Our limit on the radio afterglow from EVLA, 
combined with the inference that v c < v x . leads to a robust 
range on the circumburst density of ^0.01-0.1 cm -3 . The jet 
break time of ~2 days leads to an opening angle of 3°-8°, 
depending on the redshift and equipartition fractions, which 
leads to beaming-corrected energies of E y ~ (2-3) x 10 48 erg, 
E k « (0.3-2) x 10 49 erg, and E tot & (0.5-2) x 10 49 erg. 
This result, along with the previous jet break constraints for 
GRBs 051221A and 050724A, suggests that there may be a 
spread in true energy release, ~10 48 -10 5 ° erg for short GRBs 
(Berger et al. 2005; Grupe et al. 2006; Soderberg et al. 2006; 
Burrows et al. 2006). 

Furthermore, our optical observations provide a limit on the 
afterglow and enabled the discovery of a putative host galaxy 
with i ~ 24.3 mag. A comparison of the X-ray and optical 
data at St = 17.7 hr provides a lower limit on the host galaxy 
extinction of A 4 " st > 0.2-0. 6 mag. This is consistent with 
the high intrinsic column density from X-ray absorption when 
compared to the mean for the short GRB population. 

GRB 111020A demonstrates that rapid multi- wavelength 
follow-up is vital to our understanding of the basic properties 
of short GRBs: the geometry, energetics, and circumburst 
densities. In particular, the search for jet breaks on timescales of 
> few days is imperative for placing meaningful constraints 
on the opening angle distribution. Ideally, the detection of 
breaks in both optical and X-ray data leads to an unambiguous 
and tight constraint on the opening angle; however, optical 
afterglows are only detected in ~30% of Swift short GRBs, while 
X-ray afterglows have been detected in ~70%. Furthermore, 
optical afterglows are intrinsically faint and subject to host 
galaxy contamination, making long-term monitoring highly 
challenging. Therefore, the jet break search is optimized in the 
X-ray band, where the burst is not subject to such contamination 
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and the afterglow brightness is virtually independent of the 
typically low circumburst densities. The X-rays also allow for 
a measurement of the kinetic energy of the outflow. Deep radio 
limits provide additional constraints on the circumburst density 
and energy. The EVLA upgrade is now enabling us to probe 
events with relatively low energy scales of ~1() 48 erg and 
densities of < 1 0“ 2 cm -2 . 

The collimation of short GRBs will undoubtedly further our 
knowledge of their true energetics and rates. While the former 
provides information on the explosion and energy extraction 
mechanisms, the latter is crucial for understanding the relation 
to various progenitor systems (e.g., NS-NS mergers). Signifi- 
cant improvement on the estimated short GRB observed rate of 
~10 Gpc 3 yr" 1 (Nakar et al. 2006) will have a critical impact 
on estimates for coincident short GRB -gravitational wave de- 
tections in the era of Advanced LIGO/ VIRGO (Abadie et al. 
2010). Furthermore, a more complete knowledge of the short 
GRB redshift distribution will inform our understanding of the 
fraction of short GRBs that may originate in globular clus- 
ters, highly relevant to gravitational wave event rate estimates 
(Hopman et al. 2006; Salvaterra et al. 2008; Guetta & Stella 
2009). The uncertainty in the observed short GRB rate is domi- 
nated by the uncertainty in the beaming fraction, and with only 
a handful of short GRB opening angles measured to date, the 
discovery of even a few additional jet breaks in the coming years 
will enable significant progress. 
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